Realizing Tri-bimaximal Mixing in Minimal 
Seesaw Model with S4 Family Symmetry 

Zhen-hua Zhao 

Institute of Theoretical Physics, Chinese Academy of Sciences, 
and State Key Laboratory of Theoretical Physics, 
P.O. Box 2735, Beijing 100190, China 

^ ■ July 6, 2011 

Abstract 

\f} ' In this paper, we realize the tri-bimaximal mixing in the lepton sector 

in the context of minimal seesaw in which only two right-handed neutri- 
nos are introduced, with the discrete group 5*4 as the family symmetry. 
In order to constrain the form of superpotential, a Z3 symmetry is also 
introduced. In the model, the mass matrices for charged leptons and right- 
f> |' handed neutrinos are diagonal. The unitary matrix that diagonalizes the 

light Majorana neutrino mass matrix is exact tri-bimaximal at LO, and 
is corrected by small quantities of O(0. 01) at NLO. The mechanism to 
get the particular scalar VEV alignments used is also presented. Phe- 
nomenologically, the mass spectrum is of normal hierarchy with mi = 0, 
■ and ^2 mi and \m ee \ are about 0.058 eV and 0.003 eV respectively. 

in 



Or 



! 1 Introduction 



(~* > Nowadays, that neutrinos oscillate among three active neutrinos has been 

established as a fact by many experiments ranging from solar and atmospheric 
neutrino experiments to reactor and accelerator neutrino experiments pQ . Thanks 
to the accumulation of data, the mixing angles and mass-squared differences that 
describe neutrino oscillations now can be determined with a high precision. Ac- 
cording to a global neutrino oscillation data analysis within the three-flavor 
framework [2], the best fit values of oscillating parameters are given as fol- 
lowing, Am 2 12 = (7.59±°- 2 °) x lO~ 5 eV 2 , Amf 3 = (2.45±g-g|) x 10- 3 eV 2 and 
sin 2 (0 12 ) - 0.312±8;8& sin 2 (0 23 ) = 0.51±g;g|, sin 2 (0 13 ) = 0.010™. 

Remarkably, the mixing angles are highly consistent with that of the so- 
called tri-bimaximal mixing matrix proposed by Harrison, Perkins and Scott 
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which consists of sin 2 (6*i 2 ) = |, sin 2 (0 2 3) = \, sin 2 (0i 3 ) = 0. The particular 
oscillating pattern possessed by neutrinos might be just an accident, however, 
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if taken seriously, it maybe imply an underlying family symmetry in the lepton 
sector. Based on this point of view, many discrete groups [1], one of which is 
54 [5], have been considered as the family symmetry to study Fermion mass 
matrices. 

On the other hand, the smallness of neutrino masses is explained by the well- 
known seesaw mechanism [B] which can realize the dimension-five Weinberg's 
operator]?] by introducing heavy particles in addition to the SM particles. The 
simplest way to realize the seesaw mechanism which is compatible with observa- 
tional neutrino mass spectrum is to introduce two heavy right-handed neutrinos 
and this case is called the minimal seesaw. Ever since the celebrated paper by 
Frampton, Glashow and Yanagida [5] , minimal seesaw models assuming one or 
two zeros in the Dirac mass matrix which can lead to predictive results due to 
less parameters have been extensively studied ED ■ Interestingly, in [TT] , the 
authors make a bottom-up analysis and derive what form the Dirac mass ma- 
trix should have so that the mixing matrix is tri-bimaximal in the light neutrino 
sector after seesaw mechanism is carried out. 

In this paper, we discuss a model which can naturally give the tri-bimaximal 
mixing in the lepton sector in the context of minimal seesaw, using 54 as the 
family symmetry. In , all the special forms of the 3x2 Dirac mass matrix in 
the neutrino sector which can lead to the tri-bimaximal mixing are given, when 
the mass matrices for charged leptons and right-handed neutrinos are diagonal. 
Our model makes one form of those listed in [111 come out naturally from an 
underlying theory which has the 54 family symmetry. In the model, the mass 
spectrum is of normal hierarchy with mi =0. The mixing matrix is exact 
tri-bimaximal at LO, and is corrected by small quantities at NLO. In addition, 
the model is predictive because there is less effective parameters in the mass 
matrix than conventional seesaw model in which three right-handed neutrinos 
are introduced. 

We notice that in [T2] and [13], the authors also discuss how to get the tri- 
bimaximal mixing with 54 flavor symmetry in the minimal seesaw model. How- 
ever, their models arc different from ours in field contents, charge assignments 
and the phenomenological consequences. Besides, in [T3] the authors assume 
some particular relations among the VEVs of different Higgs fields and different 
Yukawa coupling coefficients without giving the underlying theory while in our 
model the particular forms possessed by the mass matrices come out naturally 
once suitable particle contents and charge assignments are chosen; in [13] the 
authors assume that the VEVs of scalar fields have some particular forms but 
do not discuss the scalar potential in detail while in this paper we will examine 
the scalar potential and show that the particular VEV alignment is indeed the 
minimum. 

The other parts of this paper are organized as follows: We first simply review 
some fundamental aspects of 54 which will be used throughout all the paper 
in Section 2. The model which can naturally realize the tri-bimaximal mixing 
in the lepton sector will be presented in Section 3 and some phenomenological 
consequences will also be discussed. In Section 4, we show that the particular 
VEV alignment which is essential to realizing the tri-bimaximal mixing can 
come out naturally due to the special form of the scalar potential. In Section 5, 
we discuss the impact of NLO corrections. Finally, the results are summarized 
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Table 1: Transformation properties of all the fields under 54 x Z3 x U(1)r. 



in Conclusions. 



2 S4 revisited 

54, the permutation group of 4 objects, consists of 24 elements and 5 irre- 
ducible representations which are denoted as li, 1 2 , 2, 3i, and 3 2 where the first 
number indicates dimension of the representation while the subindex is used to 
distinguish inequivalent representations of the same dimension. The presenta- 
tion of 54 we will use is the one reported in Appendix A of |14j . In the literature, 
different presentations [13 [H] have also been considered and are related to each 
other by unitary transformations which will not affect the physics. This chosen 
presentation is suitable in the sense that the mass matrices of charged leptons 
and right-handed neutrinos are both diagonal as we shall see which will bring us 
with much convenience and make the tri-bimaximal mixing more transparent. 
We only list the multiplication rules below and refer the reader to |14) for more 
details about the Clebsch-Gordan coefficients. 

with rj any representation 
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3 The structure of the model 
and phenomenological analysis 

The transformation properties of the matter fields in the lepton sector and 
fiavon fields needed are presented in Table 1. The 3 generations of left-handed 
lepton doublets are assigned to transform as 3a under 54, while e c and (/i c ,r c ) 
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transform as li and 2 respectively. The two right-handed neutrinos Ni and 
N 2 which are introduced to realize seesaw mechanism transform as li and I2 
respectively. The S*4 symmetry is then spontaneously broken by Savon fields, 
singlets under SM gauge group. The subindex of flavon fields ipi, <f>i, ip v and 
4> v indicate their roles in generating charged lepton and neutrino masses respec- 
tively. Since we discuss the model in the context of supersymmetry, two Higgs 
doublets which are singlets under S4 are included. A U(1)r symmetry is also 
introduced which will play an important role in discussing the VEV alignment 
as we shall see in Section 4. In addition, a Z3 symmetry is added which can sig- 
nificantly constrain the form of Yukawa coupling and scalar potential as shown 
in the following. As far as the quark sector is concerned, we naively assume that 
quark fields transform as singlets under S4 x Z3 and get masses via conventional 
means, so we only discuss the lepton sector in the present paper. 

The complete superpotential can be written as W= Wi + W v +Wd, 

Wi = yi e c Lh d ipi/A + y 2 {n c , r c )Lh m /A + y 3 (p c , T c )Lh d c/)i/A, (2) 

W v = vaN^K^/K + y 5 N 2 Lh u (f> u /A + M1JV1JV1 + M 2 N 2 N 2 , (3) 

where A is the cut-off scale of the theory and Wd which describes the scalar 
potential will be given in Section 5. As mentioned above, the Z3 symmetry 
constrains ipi and <j)\ to only couple with right-handed charged lcptons so that 
give the charged lepton masses, highlighting the meaning of the subindex; so 
are <p v and </>„. 

We assume that flavon fields get the following VEV configurations [T7] for 
the time being and then show it is a natural result in Section 4: 




where A, B, C, D should obey the following relations, 

9l A 2 + g 3 B 2 = 0, (6) 



C 2 = ^ M^o D 2 = 25 M^o. (7) 

2#45s + g 5 g 6 v 6545s + 3^536 

With the above specific VEV configuration, the charged lepton mass matrix 
is diagonal when the Higgs field also gets VEV: 

ViA/A 

Mi = V d I y 2 A/A + y 3 B/A ) (8) 

y 2 A/A - y 3 B/A 



where Vd is the VEV of hd- Referring to the experimental data on the charged 
lepton masses, yiVdA/A, y 2 V d A/A and y 3 V d B/A are about 0.5MeV, 900MeV 
and — 800MeV respectively. If we assume that all Yukawa coupling constants 
and all VEVs of flavon fields are both of the same order of magnitude, the hierar- 
chy between 0.5MeV and 900MeV and — 800MeV poses a problem. In order to 
explain the hierarchy, we can take into consideration the well-known Froggatt- 
Nielsen mechanism [T5] and only need to assign e c a F — N number N, so that 
electron mass becomes y±Vdt N A/ A with t, a small number, characterizing the 
ratio of the VEV of Froggatt- Nielsen flavon field to the masses of heavy Froggatt- 
Nielsen Fermion fields. From now on, we will not write Froggatt-Nielsen flavon 
field explicitly and just need to keep in mind that terms involving e c have an ad- 
ditional suppression of t , of order 10™ 3 . If we parameterize ($)/A as u where 
(<&) represents the common magnitude of flavon fields VEVs, we can get a lower 
bound for u from the requirement for Yukawa couplings to be perturbative. In 
the case of small tan/3, we can obtain u ~ 0.01 from yVd(&)/A ~ IGeV when 
y reaches order 1. 

The mass matrices for right-handed Majorana and Dirac neutrinos [pj are 



/ y 5 D/A \ 
M D = V u y 4 C/A y 5 D/A , (9) 
\ -ViC/A y 5 D/A J 

where V u is the VEV of h u . 

For convenience, we denote yiV u C ' / A and y^V u D/A as a and b. After seesaw, 
the light Majorana neutrino mass matrix is given by, 



M - ( Ml ° 1 
Mn - { M 2 J 



M v = MoM^Ml 
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(10) 



which is diagonalized by the tri-bimaximal mixing matrix with three eigen- 
values being mi=0, ^2=^- an d TO 3— in this case, m2=-\/Am^ 2 ~ 0.009 
eV, vn-i—^J Ato^ 3 ~ 0.049 eV, resulting in m, ~ 0.058 eV. Unfortunately, 
|m ee | is about 0.003eV, far below the sensitivity of present and future experi- 
ments [Til] searching for neutrino- less double beta decay. As pointed out in [TT] , 
this model cannot fulfill the role of leptogenesis [50], due to the vanishing of 
Im(MDMj-))i2. As analyzed in the charged lepton sector, yV u ($>)/A ~ IGeV in 
the case of low tan/3, so Mi and M 2 are of the order of 10 11 GeV about lower 
than in the conventional seesaw mechanism by 4 orders of magnitude. 
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4 The VEV alignment 



In this section, we show how to get the specific VEV alignment used in the 
above section. In table 1, there are several driving fields denoted by subindex 
'0' whose U(l)n charge is 2. The driving fields are gauge singlets and assumed 
to achieve no VEV. Due to the U(1)r symmetry, the superpotential terms in- 
volving driving fields has to take on the form of one driving field multiplied by 
combinations of flavon fields. In addition, the Z7> symmetry constrains <pi and 
4>i to couple with \ , and tp v and 4> v to couple with ip° and oj°. Consequently, 
the driving superpotential can be written as, 

w d = gix°(vm) + 92X°(vi4>i) + 93X°(Mi) + + 9^°(vuVu) 

+g54>°{4>»<Pv) + 96^°(<Pu<Pv) + g 7 uj Q {ip v (j) u ) + g$u>° \(j) v (j) v ). 

The conditions for the driving fields not to get VEVs can be obtained by 
deriving Wd with respect to each driving field, 



gi(2(pf 2 - 2(pnip l3 ) + g 2 (<pi 3 (j>ii ~ VuM + ff3(2^ 2 ~ 2<fe<fe) = , (11) 
gi(2<p% - 2(pnipi 2 ) + g2(cpn<f>a - ipi2<f>n) + gz{2<i>% - 2<fc<fe) = 



M^o + gi{ip 2 vl + 2(p v2 tp v z) + 55(^1 + 2<\} v2 (\) vZ ) = 

96wts + 2yV2Wl) - 97(<Pv3<f>v3 + Vul4>u2 + ip v 2<\>v\) +3s(^ 3 + 2</>„l</>„2) = .(12) 



The three equations in Eq.(ll) are satisfied by the VEV alignments for ipi 
and 4>i in Eq.(4) with A and B satisfying Eq.(6). However, the flat direction 
described by Eq.(6) poses a problem because we can only define the relation 
between A and B rather than determine their concrete values. If we include 
soft mass terms m 2 0t and m 2 ^ and assume that they are negative, (cpi) can slide 
to a large scale and then values of A and B can be fixed, as argued by Altarelli et 
al. in [IS]. The three equations in Eq.(12) are satisfied by the VEV alignments 
for ip v and <\> v in Eq.(5) and C 2 and D 2 can be determined, given in Eq.(7). 
In a word, the VEV alignment used in this paper is the minimum of the scalar 
potential. 



5 NLO corrections 



In this section, we firstly analyze the next-to-leading-order(NLO) corrections 
to Wd which lead to deviations of flavon VEVs from the LO results presented in 
Eqs.(4-5) and then discuss the NLO corrections to the mass matrices for charged 
leptons and neutrinos. Finally, we formally present how the mixing matrix is 
corrected up to the order of u. As a matter of fact, all possible £4 x Z3 invariant 
terms of dimension-six should be included at NLO. The NLO operators will be 
suppressed by a factor ($)/A ~ u with respect to the LO operators. 
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Wd can now be written as Wd = Wjj + AWd, where AW d are all possible 
dimcnsion-6 and SA x Z3 invariant terms involving driving fields. Terms con- 
tained in AWd can be obtained by inserting one flavon field to the LO driving 
superpotential and are suppressed by A. Because of the ZZ symmetry, only 
insertion of or (f> u is allowed. Hence, AWd can be written as following, 

AWd = j(£ ^f + E «if + E ^f) ( 13 ) 

where r,, Sj and U are coefficients, If , if and If arc NLO operators con- 
taining x , ip° and to respectively. The complete terms are explicitly listed 
below(we indicate the representation R with (• • •)#): 
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The corrected VEV alignment is obtained by deriving the new superpoten- 
tial with respect to the driving fields. However, the number of equations is 
smaller than the number of VEVs, we cannot solve the correction terms explic- 
itly. Nevertheless, it can be seen that the VEV alignment given by Eqs.(4-7) is 
indeed a stable minimum and the correction terms are smaller than LO terms 
by a factor u which is about 0.01 when Yukawa coupling constants are of order 
1. In the following, we will omit the corrections to VEV alignments which will 
not affect the final result. 

The NLO operators contributing to the charged lepton masses are obtained 
by inserting one ip v or <j) v to the LO Yukawa coupling superpotential and enu- 
merated in the following, 

AW i = e c Lh d [y\ 1 \^ v )+yf\^ l( t> v ) + y\ 3 \^ v )+y\ i \Mu)hJ^ 

+[( M C , T c )LM 3l [yfW,) + yfW») + 2/ ; (7 W.) + yfW^kM 2 
+ [( M C , r c )Lhdh 2 [y^ {ynp v ) + y^ + + y ; (12) (M»)hJA 2 - 

Replacing the flavon fields and Higgs field with their LO VEVs in the terms 
above, we can obtain the corrected charged lepton mass matrix. For conve- 
nience, only the order of magnitude of mass matrix entries are given in the 
following, 

/ 0(ut N ) 0{uH N ) 0(u 2 t N ) \ 
Mi = V d 0(u 2 ) O(«/10) 0{u 2 ) . (14) 
V 0(u 2 ) 0{u 2 ) 0{u) J 
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Correspondingly, the unitary matrix that diagonalizes the mass matrix MiM\ 
at NLO can be written in the form: 



/ 1 \ 

Ui = 1 V&u , (15) 
V -Vj 3 u 1 / 

where Vj 3 is coefficient of order one, V[ 2 an d V[ 3 ar e of order 0.1 and are omitted. 

Similarly, the NLO operators contributing to neutrino mass matrix are also 
obtained by inserting one ip v or <j) v and given in the following, 

AW V = Lh u N\[y^ (cp^ip^+yl 2 ^ ((p^ip^+y^ [cj) v (f> v )]/ A 2 +Lh u N 2 [y^\y> v (j) v )\/ A 2 

The modified Dirac neutrino mass matrix can be obtained by replacing the 
scalar fields with their VEVs, 

(lyP&^yPCD \ 

AM D = V u /A 2 2yl 1) C 2 +y ( l ? ) CD -3yi 4) CD . (16) 

V 2y^C 2 +yPCD Syl^CD J 
For convenience, AMp can be rewritten as, 

Ic \ 

AM D = u\ d -e , (17) 
\d e ) 

where we denote C/A, (2yl 1) C - 2yi 2) D)VJA, (2yi 1} C + y {2) D)V u /A and 
3y^DV u /A as m, c, d and e respectively. 

The NLO correction to the right-handed neutrino masses is the following, 
AM N = /liViA^^/A + f 2 N 1 N 1 <f> l ,(t> v /A + f 3 N 2 N 2 ip v ip v /A 

+UN 2 N 2 ct>vfa/A + fsNiNitpvfo/A. 

The last term vanishes when the scalar fields are replaced with their VEVs in 
Eq.(5), so there is no mixing term between JVi and N 2 up to the NLO corrections; 
as a matter of fact, if we assume that Mi and M 2 are of (9(A), then terms in 
AMpf are suppressed by u 2 with respect to terms in M^. Thus, we can disregard 
the NLO corrections to the right-handed neutrino masses. 

The modified light neutrino Majorana mass matrix can be computed as 
follows, 



AM V = MoM^AMl + AM D M^Ml (18) 



Resulting from the smallness of the NLO corrections, the unitary matrix 
that diagonalizes M v has a small deviation from the tri-bimaximal mixing. As 
a consequence, we can expand U v around Utb to the leading order of u, 



U v = Utb + AU„u, 
where AU V are parameterized by three parameters V^, V" 3 and V^, 
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Finally, the mixing matrix in the lepton sector can be obtained through the 
following relation, 



As a result, 



U = u}u v . 



sin 9 



23 



(22) 



—V&u-V^u, 


(23) 


- V&u, 


(24) 




(25) 



sin 0i2 ~ 

We can see that the NLO corrections lead to deviations of order u from the 
tri-bimaximal mixing. 



6 Conclusions and discussions 



In conclusion, with S4 as the family symmetry, we manage to obtain the tri- 
bimaximal mixing in the lepton sector in the context of minimal seesaw in which 
only two right-handed neutrinos are introduced. In order to prevent unwanted 
terms, an additional Z3 symmetry is introduced which plays an important role 
in the model as we have seen. In the model, the fields are assigned to appropriate 
transformation properties under 54 x Z3. When the flavon fields get particular 
VEV configurations, the charged leptons get a diagonal mass matrix in the 
flavor basis while the mass matrix in the neutrino sector holds a special form 
which is diagonalizcd by the right tri-bimaximal matrix. The success of the 
model heavily relies on the specific transformation properties of different fields 
under family symmetry and the special forms of the VEVs of flavon fields. 
Accordingly, whether the particular VEV configuration used in the model can 
be got naturally is a key question. Fortunately, the VEV configuration proves 
to be the minimum of scalar potential as we see in Section 4. NLO corrections 
have also been considered, and make the mixing matrix deviate from the tri- 
bimaximal pattern by quantities proportional to u whose value is about 0.01. 
Phenomenologically, the mass spectrum is of normal hierarchy with mi = 0, and 
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J2 m i an d \ m ee\ are predicted to be about 0.058 eV and 0.003 eV respectively. 
Although the phenomenological consequences are not rich, the predictions are 
definite, distinguishing the model from others. 
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